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Abstract 



04 We scrutinize the effects of non-ideal data acquisition on the tomograms of quantum states. The presence of a weight 
(— | function, schematizing the effects of a finite window or equivalently noise, only affects the state reconstruction procedure 
by a normalization constant. The results are extended to a discrete mesh and show that quantum tomography is robust 
under incomplete and approximate knowledge of tomograms. 
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1. Introduction 

Quantum technological applications require extremely 
accurate knowledge of quantum states and of the underly- 
ing quantum dynamical processes. For the application of 
fundamental principles of quantum mechanics and quan- 
tum optics to effectively foster the real-life implementation 
of quantum technologies, accurate quantum state charac- 
terization is a crucial ingredient at the interface between 
theoretical and experimental physics. Applications cover 
wide research areas ranging from nano-science to cosmol- 
ogy- 
One of the most successful quantum state reconstruc- 
tion techniques is quantum tomography [TJ [5] with its el- 
egant experimental realizations For recent reviews, 
see [10j[TI]. A large class of quantum states, expressed in 
terms of Wigner functions, can be efficiently reconstructed 
by this method, including coherent states, Schodinger cats, 
kittles, and entangled states. In quantum optics, the state 
can be directly measured by homodyne photon detection 
P2J, providing as output result the optical tomogram. 
However, the precise characterization of quantum states 
becomes problematic, both for direct measurements and 
reconstruction procedures, when noise, imperfections and 
other practical problems deteriorate the quality or limit 
the size of the data set. This entails, from an experimental 
perspective, a precise control/manipulation of the quan- 
tum system, from the source to the careful optimization 
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of the detection apparata, and from a theoretical perspec- 
tive, the refinement of mathematical inversion techniques 
for noisy data and the extension of classical tomographic 
techniques to quantum situations [U [21 IT2rfl6] . 

In quantum applications, the most used tomographic 
methods are quantum generalizations of the Radon trans- 
form p!7rtT9] . This is true both for massive particles and 
reconstruction methods based on homodyne photon detec- 
tion. The measure of optical tomograms yields complete 
information on quantum states in terms of tomographic 
probability distributions [50]. Since, traditionally, one is 
interested in other equivalent characteristics of quantum 
states, like the Wigner or the Husimi or the diagonal co- 
herent state representation function, we shall address here 
the question of the relation between these quantum state 
characteristics and our approach. 

In practice the parameters introduced in the tomogram 
have many sources of uncertainty and an efficient tomo- 
graphic measurement of the quantum state must face three 
major problems: i) a finite window (including the effects 
of detectors, entailing coarse graining and/or binning of 
the data), ii) the presence of random errors (arising both 
from the sample and the non-ideal precision in controlling 
the quadratures), and iii) a discretized "mesh" in the data 
acquisition, that affects the generation of the quadratures 
of the tomogram. 

A realistic approach must define quantum tomograms 
through a convolution with a smooth weight function S 
that spreads the marginals and reduces to the usual "clas- 
sical" Radon transform when S becomes a delta function. 
In this case one gets a kind of "thick tomography," in anal- 
ogy with signal analysis, where a similar problem arises 
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for signal detection and also requires the use of a window 
function [5TJ [55]. Notice that in a classical description, 
where tomograms are measured by the attenuation of a 
probing beam scanning a material body, this would cor- 
respond to a finite transversal thickness of the scanning 
beam. However, this is only an analogy: in the physical 
situation we are describing there is no material body, nor 
probing beams, and the thickness is in phase space. 

Our aim is to suggest extensions of the quantum tomo- 
graphic techniques by taking into account the finiteness 
of the windows, the discreteness of data acquisition and 
the effects of noise, analogously to signal analysis. Previ- 
ous work [531 HI] focused on the definition of optimization 
strategies given a set of experimental data and nonideal 
photodetector efficiency in homodyne detection (see also 
|25j for a review). 

This paper is organized as follows. In Sec. [2] we recall 
the main concepts related to the Radon and symplectic 
transforms. Section [3] is devoted to the effects of the fi- 
nite window function. In Sec. [4] and [5] we analyze how 
the tomograms are modified by noise and a discrete mesh. 
In Sec. [6] we present some numerical results for a partic- 
ular example, corroborating our analysis. Finally, in Sec. 
[7] we sketch some conclusions and outline possible future 
research directions. 



2. Preliminaries: Radon and symplectic trans- 
forms 

We shall start by studying the role of the deformation 
associated with the window weight function in the exper- 
imental tomogram and its influence on the state recon- 
struction formula. The following framework is of general 
validity, and can be applied to massive particles as well as 
photons. However, for simplicity, let p be a given quantum 
photon state, q and p the position and momentum opera- 
tors, and ip and X the local oscillator phase and quadra- 
ture in a homodyne experimental setup. The homodyne 
tomogram is given by 



W(X, ip) = Tr {p S (X — q cos <p> — psin ip)} , 



(1) 



where Tr{^4} denotes the trace of the operator A. Ex- 
pressed in terms of the Wigner function 



W(p,q) = - I (q-fflq + Z) 



d£, 



(2) 



the tomogram ([I]) reads 



dpdg W(jp, q) 
x S (X — q cos ip — p sin ip) . 



(3) 



The Radon transform of the Wigner function Efy has been 
generalized to the following symplectic, or M, transform 



W(p,q)6(X -qn-pv) dpdq, (4) 



with fi, v € R. Its complete equivalence with ([3]) is ex- 
pressed by the relation J57[ [55] 



W^X, r cosy, r sin p) = ^ W f^'^J > 



(5) 



valid for any r > 0. Equation ^ is an easy consequence 
of the fact that the Dirac distribution is positive homoge- 
neous of degree —1. The symplectic tomogram Q can be 
easily inverted by a Fourier transform 



W{p,q) 



dX dp dv 

(2^ 



W s (X,p,v)e^ x ' q ^ pu \ (6) 



which by ^ yields the inversion formula 

W{p,q) - *[dx[ drfdpW(X,p) 
(2vr) 2 J R J R+ J T 

where T = K/27rZ is the unit circle and R + = [0, +oo). 



3. Finite window 

A finite window function can be accounted for by re- 
placing the Dirac delta function in the definition of the 
tomogram Q with a suitable smearing window function 

Wi(X, (i, v) = f W(p, q)E(X- qf i- pv) dpdq. (8) 

JR 2 

Now, by Fourier transforming Eq. pj), one gets 

dX dp dv 



W(p,q) = Ns 



i{X-qn-pv) 



where 
TVs 



1 



S(-l) 



~(z)e iz dz. 



(9) 



(10) 



S(-l) 

In operator terms the state reconstruction is achieved by 



:A/S 



dX dp dv 
H3 (27r) 2 



Wl(X,fx,i>)e 



i(X-qil,-pv) 



(11) 



This is one of our central results: the inverse reconstruc- 
tion map is independent of the window function, the only 
relic of the latter being the normalization constant A/h, 
that can be fixed by the normalization of one single tomo- 
gram. Thus, no matter how involved is the shape of the 
window function, an exact reconstruction of the state can 
be always achieved by means of (|9| or ( 11 ). This property 
is very interesting because in practical cases the experi- 
mental window function S is unknown, but the result tells 
us that the exact state reconstruction is possible without 
any precise information about S. This result is also inde- 
pendent of the features of the initial state, that can even 
be nonstationary. Interestingly, this result can be extended 
to situations were noise is present. This will be done in 
the next section. 



4. Noise 



yields the symplectic tomograms 



We start from the following important observation: the 
effects of a spread function in the tomogram (|8| are de 
facto equivalent to the presence of noise, due to quan- 
tum efficiency and/or thermalization, in a homodyne to- 
mogram (e.g., if the quadrature is determined with a finite 
precision) . The origin of such noise is due to the fact that 
any counting homodyne statististics is just a sampling of 
probabilities. Namely, any statistical event can be only 
approximately sampled and in the course of repeated mea- 
surements, the outputs are always fluctuating. All this can 
be seen as a jitter in the detected statistics. 

The equivalence between such a jitter and a finite win- 
dow can be seen as follows: if there is some jitter X + £ 
on X, where £ is a random variable with a finite variance, 
one obtains 

W«> (X, fi, v) = Tv{p5{X + t- %i -#/)}. (12) 



By setting 



E[6 (X+0] =S(X), 



(13) 



where E denotes the expectation value over the random 
variable £, one gets 



fi, v) 



(14) 



and the reconstruction can be done exactly as above. 

The jitter provides extra noise and deforms the output of 
homodyne experiments by affecting the basic inequalities 
of quantum mechanics expressed in terms of tomograms. 
For this reason, it is of practical interest to take it into 
account when extracting information on the uncertainty 
relations from the measurements of tomograms. In prac- 
tice this can be done, for example, by using the output of 
photon homodyne measurements O [IT] , given as optical 
tomograms Q. The quadrature (co) variances are directly 
expressed in terms of simple integrals containing optical 
tomograms [201 and can be used to yield a noise limit 
contribution to the experimental accuracy of uncertainty 
relations in homodyne photon detection. 

Let us look at a simple but significant example of ex- 
perimental relevance. Consider the Wigner function of the 
first excited state of a harmonic oscillator (or a single pho- 
ton state) 



2(p 2 +q 2 ) 



W(p,q) 

Its Radon transform reads 
2X 2 



l e -(p 2 W)- 



-x z 



The introduction of a Gaussian window function 



E(X) 



V2iro~ 2 



^-X 2 /2<r 2 



(15) 



(16) 



(17) 



2(X 2 /r 2 + <7 2 {2a 2 + 1)) 
V^r(2a 2 + 1)5/2 

-X 2 /r 2 (2cr 2 + l) 



XC 



with r = \fji- 

W|(0, /*,!/) = 



v 2 . Note that 
2a 2 



x / 7 rr 2 (2 ( r 2 + l) 3 



(18) 



(19) 



Therefore, the presence of a window function/noise pro- 
vokes a reduction of "visibility" of the tomogram (that 
would vanish for a = 0). However, a perfect state recon- 
struction is obtained by Eqs. ^ or (11 1. 

This proves that an exact reconstruction can be ob- 
tained even in the presence of noise, provided one has very 
(in the limit, infinitely) accurate control over the position 
of the quadrature. This extends the central result of the 
preceding section. 

5. Discreteness of data acquisition 

So far we have assumed that the window function can be 
finite (thick tomograms) , but one has access to all possible 
tomograms in a continuous range of parameters. We have 
not discussed the robustness of tomography with respect to 
the discreteness of data acquisition. Let us assume that the 
experimental tomograms are gathered only on a sequence 
of discrete values of X and tp, namely 



1V,,.„, = IV j kT,m^) 



k 6 Z, to 6 



(20) 



where T > 0, 7L^ = Z/iVZ and, for convenience, N is 
an odd positive integer. If the tomograms have a limited 
bandwidth, then for sufficiently small values of T and TV -1 
one can exactly reconstruct the whole family of tomograms 
from the knowledge of the experimental ones. This is the 
content of the Nyquist- Shannon sampling theorem [3THI3T] . 
More precisely, consider the Fourier transform of W 



W(u,e) = die 



,-iwX 



2tt 



e- Uv W(X,<p). (21) 



If its support is compact and satisfies 

suppWC-D, D 
one has 



7T 7T 



T 

N 



-iX k ui e -iip m e 



where Xk = kT, (p m = 2irm/N, and 



Xd(loJ) 



1 if (u,l)€D 
if (u,l)<£D. 



(23) 



(24) 



Let us briefly show how (23) is obtained, by considering 



the less common situation of a function on the torus T: 



dip _ ; 
2tt 



Suppose f(£) = for £ £ J, with J = {-N/2, N/2) n N. 
For definiteness let us assume N odd. Then J = {—(N — 
l)/2, . . . , (N — l)/2}, and one can consider the periodic 
extension of / 



f N (e) = J2f(£-kN), 



ke'Z 



whose restriction to J coincides with /, namely, 

J(i) = fN(t)xj(e)- 



(26) 



(27) 



This property would not be true if / did not vanish outside 
J, and the phenomenon called aliasing would occur. One 
gets 

dtp 



2tt 



7(^E ei ' 



kNip 



By recalling Poisson's formula 

2lT ^2Hv- Vk) ■ 



E ei 



ikNip 



N 



(28) 



(29) 



kei 



with ip k = 2k?r/N, and setting k = jN + m with k,jEZ 
and m £ Zjv one gets 

2tt 



E ei 



fcez 
Therefore, 



A r 



E E S & ~ Vm ~ 2 ^ ■ 



f»w - jf E 



(30) 



(31) 



By plugging ( 31 ) into ( 27 1 one gets the angular dependence 



of ( 23 1 . The linear dependence is obtained analogously by 



replacing T with K and Z^r with Z. 



By Fourier inverting (j 2 T[) one gets 
/(¥>) 



E^M^) 

E /(^)^E ei 

E fifrrJSff 



<p-<Pm) 



(32) 



with 



(33) 



sin(iVa;) 
A^ sin x ' 

which is the extension of the sampling theorem to functions 
on the torus T and their discrete spectra. For a function 
g{X) on the line, the analogous, well-know formula reads 



g( x ) = ^2g( x k) sine ( 7T 



fcsEZ 



X — x k 

T 



(34) 



with X k = kT and sine a; = x~ x sin a;. 

Therefore, the tomograms are given, for iV odd, by the 
(generalized) Shannon- Whittakcr interpolation formula 



w(x,<p) = E E w * 

>m sine I 7r 



fcgZ m£2 



X — Xk 
T 



'■S'v| - * m 1 • v.Y = ?..V y - : T. 



(35) 



For N even the formula is the same, by replacing A~ by 
AT — 1 (since in that case the angular sampling points are 
N — 1 instead of N) . Remarkably, for sufficiently small T 
and A~ _1 the reconstruction of limited-bandwidth tomo- 
grams is faithful and there is no information loss. It is in- 
teresting to observe that the application of signal-analysis 
techniques to quantum tomography can be quite straight- 
forward. For example, from the mathematical point of 



view, the extension of Eqs. (|32j) and (35) of quantum to- 
mograms to a torus is very natural. 



By plugging ( 35 ) into ([T]) one can obtain a reconstruc- 
tion formula of the Wigner function in terms of the tomo- 



graphic samples ( 20 ) . By writing 



sine 7r 



X — Xk 



T 



T 
27 



-Tt/T 



-it IT 



.Mx-x k ) dw (36) 



and by performing the integration over X one gets a Dirac 
delta function S(u> + r) whose integral yields 



w(p, ?) = E E / d ^ Sn 



T 



cp-<p„ 



(27T) 2 Jo 



Tt/T 



drre ir ( Xk ~ q cos ip ~ p sin v ^ . 



(37) 



We get 



T 



(2tt)= 



Tt/T 



dr r e" 



a sin a + cos a — 1 
4Ta 2 
. a cos a — sin a 
+1 47V ' 



(38) 



with a — n(3/T. Since the imaginary part is odd with 
respect to the rotation ip — > <p> + ir, its integral over the 
torus vanishes and we finally obtain 



W(p,q) 



k£ZmeZ N T 

~ cos(a k (<p;q,p)) - 1 

4Ta k ((p;q,p) 2 
a-k(<p;q,p) sin(a k (ip;q,p)) 



V - 'Pr, 



(39) 



4Ta k {tp;q,p) 2 

with a k (ip; q,p) = Tr(X k — q costp— p sin <p)/T (and N odd). 



Let us now consider an experimental situation in which 
there is some uncertainty on the linear and angular posi- 
tion of the quadrature, that is 



'■+> 



(?) _ 



2tt 

iV 



(m + 4 2) ), (40) 



with k 6 Z and m e Zjj, where } and {Cm^} are 
two sequences of independent identically distributed ran- 
dom variables with finite standan 
Under hypothesis (p2|), one gets 



dom variables with finite standard deviations and . 



[E 



(41) 



A Fourier transform, followed by a Radon inversion (111 
yields on average a perfect reconstruction of the state. 
Therefore, quadrature uncertainties and unbiased noise do 
not affect the result. 



Let us prove Eq. (41). Consider a function g(X) on the 
line, with limited bandwidth, namely, g(u>) = for oj ^ I, 
where / = (—tt/T, +tt/T). One can write 



2n 



2tt 



(42) 



k(-_l 



By taking the expectation value, and by using Poisson's 
formula (p9| one gets 



E 



feez 



2vr 



fcez 



(43) 



Therefore, 



E 



rE» 



/ 2tt/c\ 

E 


e -i2irfc£/T' 







that, when multiplied by x/(^), yields 



E 



h^(cj) Xl (oj)] =iflf(w)E[l] 



(44) 



(45) 



because of the condition on the bandwidth of <?. Since 



E [1] = 1, we get the linear dependence of (41 1, the proof 
of the angular dependence being analogous. 

A comment seems in order. In our analysis we sepa- 
rately addressed the effects that arise from the existence 



of a finite window (or equivalently the presence of noise) 
and those that are a consequence of partial data acquisi- 
tion. The problems that arise in non-ideal situations are 
under control in two limiting cases: when the window size 
a and the sampling steps s(= T, 1/N) are well separated, 
i.e. a <C s or s <C cr. The first case can be analyzed 
through Eq. (39 1, the second one through Eqs. ([£])-( 10 ). It 
would be interesting to understand whether a single gen- 
eral formula exists, that accounts at the same time for 
the consequences of all these effects and from which both 
above-mentioned cases arise as suitable limits. Such a for- 
mula would also enable us to elucidate whether the sim- 
ple expression we obtained in terms of the normalization 



constant ( 10 ) can be generalized to the case of a finite 
number of measurements, or when noise is combined with 
other sources of uncertainty. This interesting aspect is left 
for future research. The following section is devoted to a 
partial numerical investigation of this problem. 



6. Robustness of tomograms 

So far, our analysis has taken into account the effects 
of noise and discreteness of data. We have seen that for 
sufficiently small T and A -1 the reconstruction of limited- 
bandwidth tomograms is faithful and there is no infor- 
mation loss [see Eq. (35) and following comments]. Also, 



quadrature uncertainties and unbiased noise do not affect 



the reconstruction [see Eq. (41) and following comments]. 
By combining all these results we now discuss the robust- 
ness of quantum tomograms under the afore-mentioned 
sources of uncertainty and limitations (discreteness of the 
sampling and finite precision in the determination of the 
quadrature, since we have already argued that a finite win- 
dow is equivalent to the presence of noise). Let us look 
again at the tomogram (16) and numerically investigate 



the effects that arise due to a discrete mesh and the pres- 



ence of noise on the Wigner function ( 15 ). We shall recon- 
struct the Wigner function by using Eq. (39 1 with T = 0.1, 
N = 5, and J2kez replaced by Yl\k\<K> w ^ n ^ ~ ^0. 

Figure [l] shows the tomographic reconstructions of the 
Wigner function (15), starting from its tomograms (16). 



The Wigner function in the noiseless case is shown in the 
upper panel and is practically indistinguishable from the 
original (15). Notice that, even though the tomogram 
is not band limited and the number of points K is fi- 
nite, the reconstruction by means of (39) is practically 



alias free. The reconstructions in the central and lower 



panel are affected by noise, as in (40) with zero means 
and ct^ 1 ' = = 0.5: these are large values, of the same 
order of the sampling periods T and 2?t/N (one should 
notice, however, that in this particular case the noise on 
the angular position of the quadrature does not affect the 
procedure due to the symmetry of the state considered). 
The central panel refers to a single realization of the noise, 
while the lower panel to an average over 10 realizations of 
the noise. 



5 




Figure 1: Reconstruction starting from the tomograms l |16| of the 
Wigner function W(p, q) in Eq. Jl5| . We set N = 5,T = 0.1 and K = 
40 in Eqs. ( |39[ ) and Upper panel: no noise. The reconstructed 

and the original function are indistinguishable. Central panel: single 
noise realization, cr^ 1 ' = <r^ 2 ' = 0.5. Lower panel: average over 10 
realizations of the noise. 



The reconstruction error is measured by the distance 

e = \\\W-W Tec \\ B 

= lJ B \W(p,q)-W Iec {p,q)\dpdq, (46) 

where B is the box considered (in our case B = 
[-1.5, 1.5] 2 ). We obtain e = 1.4 x 10" 4 , 0.11, 0.028 for 
the upper, central and lower panel, respectively. This 
should be compared with ||W||b — 1-01. It is apparent 
both from Fig. [l] and the above numerical results that the 
reconstruction is robust at the same time against noise and 
discretization effects. 

7. Conclusions 

In conclusion, we have discussed the practical problems 
that arise in quantum (homodyne) tomography. First of 
all, the presence of a weight function has been shown 
to introduce only a normalization constant, which entails 
no loss of information and permits the exact reconstruc- 
tion without any precise knowledge of the window func- 
tion of the experimental setup: on average, quadrature 
uncertainties and unbiased noise do not affect the recon- 
struction. Second, the discretization procedure affects the 
global reconstruction in the same way as it does in the 
classical Nyquist-Shannon setting: if the bandwidth of the 
tomogram is limited, there is no information loss for a 
sufficiently dense sampling. We have also discussed the 
most general case, in which both problems arise simulta- 
neously. Although we were not able to derive a general 
formula, from which both sources of non-ideal data ac- 
quisitions are present at the same time, and from which 
both situations investigated here are derived as limiting 
cases, we proved by numerical methods the robustness of 
the overall reconstruction against the different sources of 
imperfections. The validity of the approach proposed here 
should be tested in conjunction with other refined theo- 
retical tools based on the maximum likelihood estimation 

dum [23]. 

Since nowadays quantum states reconstruction is based 
on the measurement of tomograms, the robustness against 
non-ideal data acquisition provides solid foundation to the 
compatibility of tomographic experiments performed in 
different laboratories by different methods. Our results 
provide therefore a solid basis for the theoretical and phe- 
nomenological analysis of real tomograms. 
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